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This Letter studies the quasispecies dynamics of a population capable of genetic repair evolving 
on a time-dependent fitness landscape. We develop a model that considers an asexual population 
of single-stranded, conservatively replicating genomes, whose only source of genetic variation is due 
to copying errors during replication. We consider a time-dependent, single-fitness-peak landscape 
where the master sequence changes by a single point mutation every time r. We are able to 
analytically solve for the evolutionary dynamics of the population in the point-mutation limit. In 
particular, our model provides an analytical expression for the fraction of mutators in the dynamic 
fitness landscape that agrees well with results from stochastic simulations. 
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Genetic repair plays an important role in maintain- 
ing the integrity of nearly all organismal genomes [1]. 
Two important examples of repair mechanisms are DNA 
proofreading and mismatch repair (MMR). DNA proof- 
reading is an error-correction mechanism that is built into 
the DNA replicases themselves: If, during replication, the 
DNA replicase mistakenly pairs the wrong daughter base 
with the parent base, the replicase excises the daughter 
base and attaches a new daughter base. DNA proofread- 
ing reduces the error rate from one mismatch every 10 
bases to one mismatch every 10 7 bases [I]. 

Mismatch repair occurs immediately following daugh- 
ter strand synthesis, and requires an ability to distinguish 
between the newly synthesized daughter strand and the 
parent strand. In prokaryotes, this is achieved by exploit- 
ing the time lag between daughter strand synthesis and 
methylation of the daughter strand, while in eukaryotes 
parent-daughter strand discrimination is achieved via the 
presence of nicks in the daughter strand [TJ . 

Mismatch repair reduces the error rate by a factor of 
100, bringing the error rate in the wild- type down to 
one mismatch every 10 9 bases [TJ. Mismatch repair is 
of particular interest to evolutionary biologists, because 
it is believed that mismatch-repair-deficient strains, or 
mutators, play a crucial role in the emergence of antibi- 
otic drug resistance, and also act as gateway cells for the 
emergence of cancer [5] (and references therein). 

In static environments, mutators are generally a rela- 
tively small fraction of the population [3 , because their 
higher than wild-type mutation rate means that they ac- 
cumulate deleterious mutations more quickly than non- 
mutators. However, in dynamic environments, mutators 
can adapt more quickly than non-mutators, and so can 
rise to relatively large fractions of the population. Exper- 
imentally, it is known that in static environments the mu- 
tator fraction of Escherichia coli is on the order of 0.1%, 
while pathogenic strains under pressure to respond to an 
adaptive immune system can have mutator fractions on 
the order of 10% or more [3]. 

Because of their importance, mutators have been 
the subject of considerable experimental and theoretical 
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subject has consisted both of numerical studies modeling 
the evolution of mutation rate in time-varying environ- 
ments, as well as analytical work modeling how mutators 
influence adaptation to a new environment. 

In this Letter, we develop an analytically solvable 
model describing the influence of mutators on adaptation 
in time- varying environments. Such a model is an exten- 
sion of previous analytical work on the subject, which 
either focused on adaptation to a new, but static, envi- 
ronment [SI H]) or on adaptation in environments that 
alternate between two states [ID]. This Letter is also an 
extension of previous numerical work, in that we pro- 
vide approximate analytical expressions for the fraction 
of mutators in dynamic environments. The work here 
therefore provides a starting point for understanding the 
role of mutators in environments where co-evolutionary 
dynamics is important, such as is the case with an adap- 
tive immune system. 

An analytically solvable model describing quasispecies 
evolution in a dynamic fitness landscape was first devel- 
oped in [11] . Following the approach in [11] . we assume 
a single-peak fitness landscape in which there is one high 
fitness sequence, the master sequence a V i a ,o, in an oth- 
erwise flat landscape. We also assume that there exists 
a sequence a rePi o which corresponds to a working MMR 
mechanism. Our organisms then have genomes that may 
be denoted by a = u via a rep , where a via is the genome 
region that controls viability and cr rep is a region that 
controls repair. a V i a consists of L V i a bases, and a rep con- 
sists of L rep bases, so that the total genome length is 
L = L v i a + L rep . Note that we are essentially consider- 
ing a two-locus model, which is the approach taken by 
previous authors [SI [5J. The difference here is that we 
are assuming master-sequence-based fitness and "repair" 
landscapes, so that, in contrast to previous work, various 
transition probabilities do not need to be considered as 
independent variables, but rather may be computed from 
per-base error probabilities and sequence lengths. 

A viable organism, i.e. an organism for which o~ V i a — 
&via,0i nas a first-order growth rate constant n a — k > 1, 



2 



while for an unviable organism, k„ = 1. We define a 
repair landscape in an analogous manner: e a = e r e for 
organisms with a rep — a r ep.o and e a — e otherwise, where 
e is the per-base replication error probability and e r is the 
per-base repair error probability. 

To create a dynamic landscape, we move the fitness 
peak in the genotype space to one of its nearest neigh- 
bors (chosen randomly) at regular time intervals r. We 
then define the parameters noo, n oi, n w and riu as fol- 
lows: uqq is the number of viable non-mutators, noi is 
the number of viable mutators, nio is the number of non- 
mutators that will be viable after the next peak shift, and 
nn is t ne number of mutators that will be viable after 
the next peak shift. 

We assume that e is sufficiently small that only point 
mutations are important, and we also assume that n w 
and nn are zero immediately after a peak shift. That is, 
when the fitness peak shifts to a new master sequence, 
the number of organisms that will be viable after the next 
peak shift is assumed to be negligible. 

Based on these assumptions, we have that, during the 
time interval [nr, (n+l)r) immediately following the n^ n 
peak shift, the quasispecies equations take the form, 
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We may solve Eq. (1) for t 6 [nr, (n + l)t) with the 
initial conditions nio(nr) = nn(nr) = 0. Noting that at 
t = (n + l)r the fitness peak shifts to the next master 
sequence, we have, 



nao((n + l)r) = A w (T)n ao (nT) 

n i((n + l)r) = Ai 1 (r)n 00 (nr) + Su(r)n i(nr) 
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where the coefficients Aiq(t), Ah(t), Bu(t) are given 
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Let us denote by r n the ratio of viable mutators to 

th 

viable non-mutators immediately following the n peak 
shift. Then from the solution given above we have, 



' n+l 



a(r) + /3(r)r v 



(•5) 



where a(r) = A x1 {t) / A w {t) and (3(t) = B u (t)/A w (t), 
from which it can be shown that r n = a(r) ^2™ = q [P(t)]' 1 + 
f3 n ( T ) r o- If < /? < 1, then, as n — > oo, we obtain 
that ee lim n ^ oc r„ = AioCt^-bLct) ' wmch implies a 
periodic solution where the fraction of viable organisms 
that are non-mutators at the beginning of each cycle is 
nonzero. When [3 > 1 we obtain that r n diverges, and so 
the fraction of viable organisms that are mutators goes to 
1. This phenomenon is known as the repair catastrophe 

For the periodic solution that develops after a suffi- 
cient number of iterations, the fraction of viable organ- 
isms that are mutators at the beginning of a peak shift is 
then r ^/(l + r^x,). As r — * 00, it may be readily shown 
that the fraction of viable organisms that are mutators 
just before the next peak shift is given by the static land- 
scape expression obtained in jl]. 

In order for the population to form a stable quasis- 
pecies and remain viable in the time-varying landscape, 
the relative growth of the viable population between the 
peak shifts should be larger than that of the background 
(genome sequences some distance away from the master 
sequence). Since the growth constant of non- viable or- 
ganisms is equal to 1, the criterion for adaptability is 
given by, 

n o((n + 1)t) + n ai ((n + l)r) ^ ^ T 
noo(nT) + noi(nr) 

If this condition is not met there will be a drift from the 
master sequence to the "background" , corresponding to 
the error catastrophe. 

When the fraction of viable organisms that are non- 
mutators is nonzero, Eq. (6) may be shown to be equiv- 
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while when the fraction of viable organisms that are non- 
mutators is zero, i.e. when the population has under- 
gone the repair catastrophe, Eq. (6) may be shown to be 
equivalent to, 
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We may define three distinct parameter regimes, each 
corresponding to the viability or non- viability of the non- 
mutator and mutator populations, respectively. 

The first region is defined by K nm > 1, n m . Here, the 
non-mutators' effective growth rate is both larger than 
that of the background as well as that of the mutators, 
so that a stable quasispecies is formed with a non-zero 
fraction of viable non-mutators. 

The second region is defined by n m > 1 , K nm . Here 
the mutators' effective growth is both larger than that of 
the background as well as that of the non-mutators, so 
that a stable quasispecies is formed, but it consists en- 
tirely of mutators. The transition between the first and 
second regions corresponds to a localization to dereal- 
ization transition over the repair portion of the genome 
that is termed the repair catastrophe [2J. 

Finally, the third region is defined by n nm , n m < 1, so 
that no stable quasispecies forms, and the population is 
unable to adapt to the changing fitness landscape. 

The solutions to the equations K nm ,K m — 1 set the 
upper and lower mutation thresholds for the population 
to survive (see Figure 1). The upper mutation threshold 
is the ordinary error catastrophe for the static landscape, 
while the lower threshold appears only in the case of a 
dynamic landscape. This lower threshold arises because, 
in a dynamic fitness landscape, the mutation rate must 
have a minimal value to allow the population to adapt. 
Below this dynamic error threshold, the population is 
unable to adapt, and so no stable quasispecies is formed 

nu. 

There is some minimal value of r, denoted T min , below 
which the population cannot adapt. This corresponds to 
a threshold rate of change of the fitness landscape above 
which the population cannot adapt, irrespective of the 
error threshold. This minimal value for r arises because, 
as t is decreased, the mutation rate must increase to 
maintain a stable quasispecies. Once this minimal mu- 
tation rate exceeds the static error threshold, then no 
quasispecies can exist. 

An analytical approximation for r m i n can be found by 
expanding Eq. (8) to second-order in e, solving for e, and 
then solving for the value of r where the two solutions 
are equal. Assuming that k is sufficiently larger than 1, 
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Knm = 1 and n m = 1 as a function of cycle length r. The 
area confined by the symbols is the survivability region of the 
population. In the region above the solid curve the fraction 
of non-mutators is nonzero, while below the solid curve only 
mutators exist. r min found from Eq. (9). 
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and that r is sufficiently large, so that e 
obtain the following approximate expression for t„ 
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An e versus r plot of the solution to the equation 
Knm = Km yields a continuous curve, above which it can 
be shown that n nm > n m , and below which it can be 
shown that K nm < K m . Therefore, the region above this 
curve that also lies between the lower and upper mu- 
tational thresholds of the non-mutators corresponds to 
a parameter regime where there is a stable quasispecies 
with a nonzero fraction of non-mutators. The region be- 
low this curve that also lies between the lower and up- 
per mutational thresholds of the mutators corresponds to 
a parameter regime where there is a stable quasispecies 
consisting entirely of mutators. 

We developed a stochastic code to simulate the dynam- 
ics of a population of self-replicating genomes of the form 
c = (Tvia&rep- At each time step, a given genome repli- 
cates with probability K a At, where K a is the growth con- 
stant and At is the size of the time step. We choose At 
to be sufficiently small so that the replication probability 
is much smaller than one. After every time interval r, a 
shift occurs and a new master sequence is generated to be 
one point mutation away from the old one. The mutation 
is generated at a different random location in the viable 
region to avoid population accumulation (to meet the ini- 
tial conditions of the analytical model), i.e. the mutation 
position repeats itself every L via peak shifts. Each real- 
ization simulated at least 10 4 time steps to make sure a 
periodic solution was reached. 
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FIG. 2: N = 1.5 x 10 5 , L via = 20, L r 



10, S = 4, e = 0.02, 



e r = 0.1, r = 2: A single trajectory from a stochastic simula- 
tion. Note that the mean fitness reaches its static landscape 
value. The simulated fraction of mutators at the beginning 
and at the end of each cycle shows good agreement with the 
analytical results. 



"-oo - * "oi, "oo - * "id an d "oi - * "n transition proba- 
bilities are all first-order in e, while the transition prob- 
ability n 00 — > "n is second-order in e, and is therefore 
neglected in a model that only considers mutation prob- 
abilities that are up to first-order in e. 

For larger e, the non-mutator to mutator transition 
probability, which is 1 — (1 — e r e) Lrep , goes from being 
first-order in e to zeroth-order, so that the uqq — * "n 
becomes first-order in e as well. Therefore, if we wish 
to use a first-order model for larger e, then this addi- 
tonal transition probability must be included. Including 
the noo — > "n transition for larger e significantly im- 
proves the theoretical prediction of the mutator fraction. 
For future research, we wish to develop a more system- 
atic approach for estimating the mutator fraction in a 
dynamic fitness landscape that does not require the as- 
sumption of small e. Furthermore, we would like to move 
beyond mean-field descriptions of the evolutionary dy- 
namics and develop analytical approaches for quantifying 
various stochastic effects. 
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FIG. 3: N = 1.5 x 10 5 , L via = 20, L rep = 10, S = 4, e = 0.02, 
e r = 0.1: Fraction of mutators as a function of r. The results 
are averaged over 10 independent runs. 
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As can be seen from Figures 2 and 3, as long as e is 
sufficiently small, we obtain good agreement between the- 
ory and simulation. However, due to stochastic effects, 
the fraction of mutators requires a larger population size 
than other evolution parameters (e.g. mean fitness, mean 
Hamming distance) before reasonable agreement with the 
theoretical results is obtained. 

For sufficiently large e, the agreement between our 
point-mutation model and the simulation results breaks 
down. The reason for this is that, when e is small, the 



